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CHAPTER  1 
INTRODUCTION 

Future  communications  systems  will  require  informa¬ 
tion-carrying  capacities  of  several  hundred  Mbits/sec  over 
long  link  lengths.  Graded-index  fibers  have  assumed  consid- 
""  erable  importance  since  they  offer  the  possibility  of  multi- 

mode  propagation  in  a  relatively  large  core  fiber  with  low 
mode  dispersion. 

> 

The  performance  of  optical  fibers  for  communications 

» 

is  chiefly  determined  by  optical  loss  and  by  the  temporal 
dispersion  of  light  pulses  caused  by  material  and  modal 
dispersion.  Material  dispersion  occurs  because  the  index 
of  refraction  varies  as  a  function  of  wavelength.  Modal 
dispersion  results  from  the  differences  in  the  path  lengths 
of  the  various  transmission  modes.  Thus,  material  and  modal 
dispersion  determine  the  bandwidth  and  the  information  carry¬ 
ing  capacity  of  the  cable.  For  single  mode  step-index  fibers, 
the  diameter  of  the  fiber  is  small,  which  causes  difficulty 
in  splicing  and  coupling  to  light  sources.  For  a  multi- 
mode  step-index  fiber  with  a  large  diameter  core,  modal 
dispersion  severely  limits  the  bandwidth  and  the  fiber's 
utilization  over  long  distances  [1]  . 


By  varying  the  index  of  refraction  across  the  radius 
of  the  fiber  core,  it  is  possible  to  compensate  for  the  path 
length  differences  of  the  many  modes  carrying  power  in  the 
fiber.  Graded-index  fibers  thus  combine  the  advantages  of 
large  core  diameters  with  low  modal  dispersion. 

In  this  report,  we  will  examine  the  effects  of  pulse 
dispersion  and  noise  in  a  graded-index  fiber  optic  communi¬ 
cation  system  which  utilizes  pulse  amplitude  modulation. 

It  will  be  shown  that  pulse  dispersion  can  result  in  inter¬ 
symbol  interference  and  degraded  error  performance  of  the 
system.  Expressions  for  rms  output  pulse  width,  intersymbol 
interference,  optimum  receiver  threshold  levels,  and  error 
probabilities  will  be  developed  in  terms  of  system  parameters. 
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CHAPTER  2 

PULSE  DISPERSION  IN  GRADED- INDEX  FIBERS 

2 . 1  WKB  Solution  for  Graded-Index  Fibers 

The  propagation  constants  of  a  cylindrically  symmetric 
waveguide  with  a  arbitrary  index  profile  can  be  obtained 
using  the  Wentzel-Kramers-Brillouin  (WKB)  approximation  which 
is  widely  used  in  quantum  mechanics  [2] .  Of  primary  interest 
will  be  the  rms  output  pulse  width  which  will  be  a  function 
of  both  the  material  and  modal  disperson  in  the  optical 
fiber.  In  the  WKB  analysis,  the  following  assumptions  are 
made : 

1.  The  fiber  core  diameter  is  large  so  that  numerous 
modes  will  be  excited. 

2.  The  index  variation  is  small  over  distances  of  a 
wavelength. 

3.  All  propagating  modes  carry  equal  amounts  of 
energy . 

The  refractive  index  profile  of  the  fiber  can  be 
modeled  by 


"(r)  =  nQ 

1  - 

\a)  j 

r  <  a 

(2.1) 

n(r>  =  no 

[1  - 

2A]1/2 

r  >  a 

(2.2) 

where 


r  =  distance  from  the  axis  of  the  fiber 

a  =  radius  of  the  core 

2  2  2 
A  "  (no  -  nl>/2ni 

nQ  =  refractive  index  at  the  fiber  axis 
n  •]_  =  refractive  index  of  the  cladding 
a  defines  the  index  gradient  of  the  fiber. 

It  has  been  shown  [3]  for  the  above  case  that  the  total 
number  of  guided  modes  is  given  by 


where  k  is  the  optical  wavenumber, 
a  given  mode,  m,  is  given  by 


1  (2.3) 

The  group  velocity  for 


(2.4) 


where  3m  is  the  propagation  constant  of  the  m'th  mode  in  the 
direction  parallel  to  the  fiber  axis.  The  time  taken  for 
energy  propagating  in  the  m'th  mode  to  travel  a  distance  L 
is  therefore 
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By  using  the  WKB  approximation  it  can  be  shown  that  [3] 


S>m  =  n  k 
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Therefore,  the  propagation  delay  can  be  written 


T  = 
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where 
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and  the  group  index  is  defined  as 


N  =  n  +  A  st-t— 

o  IdX 


(2.8) 


(2.9) 


X  =  wavelength  of  light  transmitted. 

The  mode  number,  m,  used  above  is  related  to  the  double 
quantum-number  LP„,  ,  system  of  labeling  by  the  simple  approx¬ 


imation 


m  -  ( 2m '  +  S, '  +  1 ) 


(2.10) 


For  the  purposes  of  specifying  the  fiber  bandwidth 
for  a  digital  system,  the  rms  pulse  width  is  of  primary 
importance  and  is  given  by 


T  2  2 

o  -  <  T  >  —  <T> 

!  m  m 


(2.11) 
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(2.12) 


where  is  the  term  due  to  intermodal  dispersion  and  is 
given  by 
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=  Si  ~  2  -  X 
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(2.14) 


3a  -  2  -  2y 
2  2 (a+2) 


(2.15) 


The  a..  term  is  due  to  intramodal  dispersion  and  equals 


where  the  rms  spectral  width  of  the  source  is  defined  by 


dA  (A- 


1/2 


'V 


S  (  A) 


(2.17) 


and  the  mean  source  wavelength  is  given  by 

00 

A Q  =  J  d  A  A  S ( A )  .  (2.18) 

o 

The  rms  pulse  width  has  been  separated  into  an  inter- 
modal  and  an  intramodal  component.  The  intermodal  term 
results  from  delay  differences  among  the  various  modes  pro¬ 
pagating  in  the  fiber.  The  intramodal  component  represents 
an  average  of  the  pulse  broadening  within  each  mode  and 
arises  from  material  and  waveguide  dispersion.  With  the 


I 
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proper  choice  of  i  in  a  graded-index  fiber,  the  intermodal 
dispersion  is  small  and  the  intramodal  effects  caused  by 
material  dispersion  predominate.  The  lower  limit  of  the  out¬ 
put  rms  pulse  width  is 


°min  — 


(2.19) 


2 . 2  Selection  of  the  Optimum  Gradient  Index 

The  total  rms  pulse  spreading  can  be  calculated  [2] 
using  (2.12),  (2.13),  and  (2.14).  In  Figure  1,  this  rms 


Figure  1.  Calculated  rms  pulse  spreading  in  a  graded  fiber 
as  a  function  of  the  parameter 
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pulse  width  is  shown  as  a  function  of  a  for  three  types  of 
GaAs  sources  all  operating  at  900  nm  but  having  different 
spectral  bandwidths.  The  fiber  parameters  are  characteris¬ 
tic  of  TiC>2  doped  silica  guide  with  a  =  0.0059.  The  three 
curves  correspond  to  an  LED,  an  injection  laser,  and  a  dis¬ 
tributed  feedback  laser  having  typical  ritis  spectral  widths 
of  15  nm,  1.0  nm,  and  0.2  nm  respectively.  For  comparison, 
the  dashed  curve  represents  the  predicted  rms  width  if 
material  dispersion  and  intramodal  broadening  are  ignored. 

In  this  case,  (2.13)  alone  is  used  to  compute  pulse  spread¬ 
ing.  The  figure  demonstrates  that  pulse  dispersion,  with 
an  LED  source,  has  a  lower  limit  of  approximately  1.0  ns/km 
and  is  set  by  the  material  dispersion  of  the  fiber.  Using  a 
laser  source,  the  pulse  dispersion  approaches  0.02  ns/km. 

In  order  to  minimize  the  pulse  dispersion,  an  optimum 
index  gradient  must  be  selected  which  will  decrease  the 
intermodal  dispersion  to  a  minimum.  In  (2.13),  if  we  set 
y=0,  which  corresponds  to  ignoring  material  dispersion,  the 
minimum  dispersion  occurs  for  ol  of  nearly  2.  For  y^O.O,  the 
minimum  occurs  at  approximately  2.23  and  is  given  by 


a  .  =  2  +  y  - 

mm  J 


/  (  4-t-y )  (3+y)  \ 
5  +  2y  j 


(2.20) 


Since  the  optimal  a  varies  strongly  with  the  wavelength  of 
the  source,  the  output  RMS  pulse  width  is  dependent  upon 
the  source  wavelength  for  a  particular  cable. 
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2 . 3  Accuracy  of  the  WKB  Analysis 

Although  the  WK3  analysis  is  relatively  simple  and 
comprehensive,  it  is  an  approximation  and  results  in  errors 
for  modes  close  to  cutoff  [4] .  In  reality,  these  modes 
always  carry  less  power  than  the  modes  that  are  well  guided. 
This  fact  means  that  predicted  values  of  rms  pulse  width 
will  not  always  correspond  to  experimental  values.  In  the 
WKB  analysis,  the  maximum  delay  difference  between  propagat- 
ing  modes  was  approximately  TA  /2,  where  T  and  A  denote  the 
propagation  time  and  the  relative  refractive  index  respec¬ 
tively.  By  more  exact  analysis,  we  find  that  the  delay  time 
differences  between  well-confined  modes  and  modes  close  to 
cutoff  approach  TA.  Since  we  have  assumed  a  small  A,  in 
this  case  0.0059,  inaccuracies  of  this  magnitude  can  be 
ignored  for  cable  lengths  under  discussion. 

In  addition,  the  WKB  analysis  predicts  an  optimal 
index  gradient.  In  practice,  it  is  difficult  to  produce 
an  optimal  index  gradient  in  a  fiber  many  kilometers  long. 
For  the  LED  source,  pulse  broadening  of  less  than  1.5  nsec/ 
km  can  be  achieved  if  a  is  within  25%  of  the  optimal  valra. 
For  the  injection  laser,  an  2  within  5%  of  the  optimal 
value  will  give  widths  less  than  0.2  nsec/km,  and  for  the 
distributed  feedback  laser,  widths  of  0.05  nsec/km  are  pre¬ 
dicted  if  1%  control  on  a  can  be  achieved. 
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2 . 4  RMS  Output  Pulse  Width  as  a  Function  of  the  Fiber 
Parameters  — 

The  temporal  spreading  experienced  by  a  pulse  of  light 
as  it  propagates  along  a  graded  index  optical  fiber  is 
characterized  in  terms  of  the  dispersion  parameters  and 
given  in  (2.13)  and  (2.16).  We  can  express  the  output 
pulse  width,  T  ,  •'■n  terms  of  the  individual  fiber  para¬ 
meters  as 

1/2 

(Tout’rms  '  [<zt,I|2  +  <ZV2  +  Ti  ]  <2-21 

where  z  is  the  fiber  length  and  is  the  rms  input  pulse 
width.  With  the  optimal  choice  of  a  for  the  graded  index 
fiber,  the  intermodal  dispersion  is  small  and  the  dispersion 
term,  predominates.  Therefore  the  output  pulse  width 

can  be  expressed  as 


(T  ) 
out  rms 


i  12  V 2 


(2.22) 


The  expression  for  o  ,  (2.16),  can  be  evaluated  using  data 
listed  in  [2] . 


In  Figures  2  and  3,  the  rms  output  pulse  width  is  plot¬ 
ted  as  a  function  of  the  fiber  length,  z,  for  various  values 
of  rms  input  pulse  width.  A  GaAs  LED  source  with  a  line 

width  of  15  nm  and  wavelength  of  900  nm  is  assumed.  When  T^ 

2 

is  significantly  larger  than  20  nsec,  the  T^  term  is  dominant 
in  (2.22)  and  TQut  increases  very  slowly  with  increasing  z. 


100  ns 


(su)  ‘hiqim  3S"ind  indino  sm 


•t4 

0) 

44 

<u 

o 

3 

c 

*■4 

3 

(0 

44 

> 

nj 

0) 

3 

t n 

0 

ns 

•H 

. — . 

3 

cn 

> 

■3 

• 

C 

S4 

, — . 

0 

0 

m 

U 

44 

•a 

a) 

c 

in 

0 

0 

in 

u 

c 

u 

(U 

<TJ 

0) 

in 

C 

44 

0 

a) 

c 

s 

nj 

■P 

0 

c 

3 

• — - 

0 

■H 

* 

■H 

' — ' 

•>. 

G 

N 

■4 

44 

G 

T3 

44 

•H 

G 

r0 

? 

44 

**4 

Cn 

s 

a> 

G 

w 

<D 

<u 

»-4 

in 

3 

<-4 

04 

>4 

3 

a) 

04 

JJ 

G 

3 

•>4 

44 

04 

44 

3 

4J 

04 

3 

aj 

G 

0 

r* 

*H 

44 

CA 

cn 

s 

44 

e 

ns 

0 

S4 

s 

S3 


In  Figures  4  and  5,  TQut  is  again  plotted  as  a  func¬ 
tion  of  2  and  TV.  A  GaAs  laser  source  with  a  wavelength 
of  900  nm  and  a  line  width  of  1.0  nm  is  assumed.  The  temporal 
spreading  experienced  as  a  function  of  increasing  z  is 
similar  to  che  results  obtained  with  the  GaAs  LED  source. 

The  magnitudes  of  TQut  and  T^  differ  in  these  two  cases 
because  of  the  larger  source  line  width  of  the  LED  source. 


CHAPTER  3 


ANALYSIS  OF  AN  OPTICAL  DIGITAL  RECEIVER 

3 . 1  Mathematical  Model  of  the  Receiver 

In  fiber  optical  systems,  the  detector  output  pro¬ 
duced  by  the  optical  power  impinging  upon  an  avalanche 
photodiode  is  an  inhomogeneous,  non-stationary  shot  noise 
process  [5].  It  is  assumed  that  the  optical  power,  p(t), 
impinging  upon  -he  av..lanche  photodiode's  p-n  junction  is 
a  truncated  street  of  independent  binary  digits.  A  mathe¬ 
matical  model  of  an  optical  digital  receiver  is  depicted  in 
Figure  6.  St  consists  of  a  p-n  junction  followed  by  a  ran¬ 
dom  multiplier  with  a  known  electron  output  probability 
density  function  followed  by  a  filter.  Arrival  or  absence 
of  a  transmitted  pulse  is  decided  by  a  comparator.  An  in¬ 
correct  decision  may  be  caused  by  thermal  noise,  intersymbol 
interference,  shot  noise,  or  dark  current. 

The  current  x(t)  is  comprised  of  xg(t),  the  signal 
current  from  the  random  multiplier  output,  and  the  Gaussian 
noise  current  xfc^(t).  x(t)  is  given  by 

x(t)  =  xs<t)  +  xth(t)  *  (3.1) 

xs(t)  is  an  inhomogenous ,  non-stationary  shot  noise  process 
given  by 


Lgure  6.  Mathematical  model  of  an  optical  digital  receiver. 
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V(t) 

Xs(t)  =  ^  gv  3(t-tv)  (3.2) 

V=1 

where  t  is  the  time  release  of  the  vth  electron  from  the 
p-n  junction.  V(t)  is  a  random  number  of  these  electrons 
which  occur  in  a  time  interval  up  to  t  and  are  statis¬ 
tically  independent  multiplicative  random  variables  correspond¬ 
ing  to  the  random  multiplication  and  random  gains  encountered 
in  the  multiplier. 

The  p-n  junction  produces  primary  electrons  in  response 
to  the  impinging  optical  power,  p(t),  and  due  to  thermal 
effects.  The  rate  at  which  primary  electrons  are  produced 
obeys  Poisson  Statistics  with  an  average  rate  of 

A(t)  =  As(t)  +  Ad  (3.3) 

where 

Ag (t )  =  ^  p(t)  (3.4) 

and 

n  =  junction  quantum  efficiency 
h  =  Planck's  constant 
f  =  optical  frequency 
>d  =  dark  current  "counts"  per  second. 

A ( t)  is  only  the  average  rate  at  which  electrons  are  pro¬ 
duced.  The  probability  that  exactly  k  electrons  are  produced 
in  any  interval  T  seconds  long  is  given  by 
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P[k, 


WV1 


(3.5) 


where 

t  +T 
o  s 

X ( t ) dt  .  (3.6) 

Given  p(t),  the  number  of  electrons  produced  in  any 
interval  is  statistically  independent  of  the  number  produced 
in  any  other  adjacent  interval  [6].  p(t)  is  assumed  to  be  of 

the  form  of  a  digital  pulse  stream  given  by 

00 

P(t)  =  d.  <  |S0(t-iTs) |2 >  (3.7) 

i  =  -00 


where 

d^  =  sequence  of  independent  equiprobable  data  symbols 
Ts  =  signaling  interval 
2 

<|S  (t-iT  ) |  >  =  fiber  distorted  Dulse  shape  which 
os' 

is  positive  for  all  t. 

We  shall  assume 

<j  S  (t-iT  )  \2>  =  1  .  (3.8) 

o  s 

Thus,  d^  is  the  energy  in  the  ith  pulse.  We  will  assume 

that  the  fiber  distorts  the  pulse  enough  so  as  to  cause 

intersymbol  interference  (ISI).  In  the  presence  of  ISI, 

X  (t)  can  be  broken  into  two  parts: 
s 
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>  I 

f 


k 


i 


L 


Ut) 


00 

-  ET  w0  <|so(t),2>  +  .X  di 


1  =  -® 
i^O 


s0(t-iTs)  r>] 


(3.9) 


+  V 


The  zeroth  term  in  the  summation  represents  the  photon 
arrival  rate  for  the  pulse  under  detection.  The  signaling 
slots  containing  the  ISI  components  are  comprised  of  the 
remaining  terms. 


3 . 2  Statistics  of  the  Receiver  Output 

In  Figure  6,  the  current  x(t),  is  linearly  filtered. 
h(t)  is  the  impulse  response  of  the  filter,  and  y(t)  is  the 
filter  output.  The  signal  component  of  this  output  (due  to 
the  nonuniform  shot  noise)  is  given  by 


y_ (t)  =  x  (t)  *  h ( t)  =  I  ah ( t-t  )  . 

5  5  *  i  n 

n=-do 


(3.10) 


The  mean  and  variance  of  the  signal  at  the  receiver  output 
can  be  calculated  using  Campbell's  theorem  [7] : 


Ely  .  ( t) )  =  A(  :)  *  h(T)  = 


) 


y  dt  Hf  p(t)  h(t" 

—  00 

Cy  (tl-c2>  ’/  d-  CT  p<Tlh(trT)h(t2-T) 

S  -TO 

where  C  is  the  covariance  of  the  received  signal.  The 
■‘s 

impulse  response  of  the  filter  is  described  by 


(3.11) 


3.12) 
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h  (t)  = 


h  i-i  <2^ 


0  elsewhere 


(3.13) 


h(t  -  t)  =  h  rect 


(3.14) 


where  h  is  a  constant .  The  region  over  which  h(t-T)  is  non¬ 
zero  satisfies  the  inequality 


1  t  -  T  ■  <  — 

I  u  1  O 


(3.15) 


which  redefines  the  limits  of  integration  to 


Ts  Ts 

~  2  *  +  ~  • 


(3.16) 


By  using  the  above  results,  the  mean  and  variance  of  y  (t) 


can  be  written  as 


E[yg(t) 


T 

^i2 


dT  KT 


(3.17) 


Var[ys(t)l  = 


dr  HT“  P(T)  =  h  E[ys(t)l 


:3.i8) 


J  T= 
‘'I2 


where  p(r)  is  given  by  (3.7).  If  the  sequence  of  data  symbols 
is  known,  the  mean  and  variance  of  the  signal  component  of  the 


output  become 
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t  +  T  /2 

r  _  M 

ECyg(t)  !di]  -  j  dr  ^  d.  <;So(-:-iTs)  ;2>  (3.19) 

j  i=oo 

-  V2 


Var [y  ( t) |d.)  =  h  E [y  (t) |d.] 

S  X  S  X 


(3.20) 


The  thermal  noise  component  of  the  output  is  y  .  (t)  where 


oo 

yth(t)  ~  f  dT  xth(T)  h ( t  -  T)  . 

—  oo 


(3.21) 


The  thermal  noise  is  uncorrelated  with  zero  mean.  The  vari¬ 
ance  of  the  filtered  thermal  noise  at  the  receiver  output  is 


"th2  -/  at  ir  h2,t> 


(3.22) 


which  reduces  to 

m 

A  s  _2 

9  r  N  9  n  t  h 

Jth  7  dt2 2  h2  =  -V—  (3-3) 

o 

where  Nq/2  is  the  two  sided  power  spectral  density  of  the 
thermal  noise. 


3 . 3  Temporal  Pu ’.se  Slot  Representation 

Our  receiver  performance  analysis  assumes  that  only  the 
pulse  slots  immediately  adjacent  to  the  pulse  slot  under 
detection  will  contain  significant  ISI  contributions. 

Figure  7  illustrates  the  temporal  pulse  slot  representation 
used  where  Tsq  contains  the  energy  of  the  pulse  ur%'  r 
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Time 


Figure  7.  Temporal  pulse  slot  representation. 


detection.  The  energy  in  Ts+1  is  integrated,  sampled  and 
dumped  exactly  one  period  before  the  energy  in  Tso  is 
analyzed.  Exactly  one  period  after  the  energy  in  Tsq  is 
integrated,  sampled  and  dumped,  the  energy  in  Ts_^  is 
analyzed.  In  Figure  7,  Iq  equals  the  energy  in  Ts  due  to 
the  pulse  under  detection.  and  I+^  are  the  energies  in 

Tsq  due  to  the  temporal  spreading  of  pulses  in  Ts_^  and  Ts+^ 
respectively.  The  1^  energies  are  described  in  integral 
form  by 


t 


! 


L 
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I  = 


L-1 


Ts 

+  2* 


T 

t  -  s 


2 

T 


t  + 


t  - 


dr  dQ  <«So(t) j 2> ! 


(3.24) 


t=0 


dT  d_x  <'So(T  +  Ts)  i*> 


(3.25) 


t=-T 


>1 


t  + 


t  - 


dT  d+1  <|SQ( 


V 


T 

_a 

2 


(3.26) 


t=T 


3 . 4  Expression  for  the  Probability  of  Error 

Amplitude  shift  keying  is  chosen  to  modulate  the  light 
intensity  in  the  fiber.  The  linear  receiver  filter  and  samp¬ 
ling  circuitry  comprise  an  integrate-and-dump  detection 
scheme  [8).  The  signal  y(t),  in  Figure  6,  is  sampled  every 
Ts  seconds,  where  Tg  is  the  temporal  pulse  slot  width,  which 
is  inversely  proportional  to  the  bit  rate.  The  sampled  sig¬ 
nal,  y(Tg) ,  is  the  sum  of  a  signal  component  s,  and  a  noise 
component  n.  s  is  due  to  the  signal  induced  shot  noise  in 
time  slot  Tsq  when  a  1  is  transmitted.  n  is  comprised  of 
thermal  noise,  and  the  shot  noise  components  spilling  over 
from  adjacent  time  slots  Ts_1  and  Ts+1-  The  dark  current 
primary  electron  count  is  assumed  to  be  small,  and  is 
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therefore  neglected  in  this  analysis.  y(T  )  is  then  compared 
to  a  decision  threshold,  r^.  If  y(Ts)  is  greater  than  rT, 
the  receiver  decides  a  1  was  transmitted.  If  y(T  )  is  less 
than  rT,  the  receiver  decides  a  0  was  transmitted. 

We  will  temporarily  neglect  ISI  effects  and  examine  the 
two  ways  in  which  errors  occur.  If  a  1  is  transmitted,  an 
error  occurs  if 

s  +  n  _<  rT  (3.27) 

where  n  is  comprised  of  shot  noise  and  thermal  noise.  An 
error  will  also  occur  if 


n  > 


(3.28) 


when  a  0  is  transmitted.  In  this  case,  n  is  composed  only 
of  the  thermal  noise  present  in  Tsq.  Since  0's  and  l’s  are 
equally  likely,  the  average  probability  of  error  can  be  given 
by 


P  (e)  =  P  (e  1 1)  •  P  ( 1)  +  P  (e  i  o)  *P  (o) 

=  j  P(s  +  n  <  rT)  P (n  >  r?) . 


(3.29) 


Our  error  analysis  concentrates  on  the  energies  in  time 
slots  Tsq,  Ts_^,  and  Ts  +  ^,  which  give  rise  to  eight  possible 
probability  of  error  combinations  based  upon  0’s  and  l’s 
appearing  independently  in  the  three  time  slots  under  dis¬ 
cussion.  The  total  average  probability  of  error  can  be 
expressed  as 
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P(e)  = 

P(e;d_1  =  1, 

o* 

+ 

>— 

ii 

H* 

•  P(d_1  =  1, 

d+l 

=  1) 

+ 

P  (e i d_^  =  0, 

d+i  =  1) 

•  P(d_x  =  0, 

d+l 

=  1) 

+ 

P(e!d_x  =  1, 

d+i  =  1) 

•  P(d_1  =  1, 

d+l 

=  0) 

+ 

P(e|d_1  =  0, 

o 

II 

H 

+ 

•  P(d_1  =  0, 

d+l 

=  0) 

Since  cL  is  a  sequence  of  independent  eauiprobable  data 
symbols  and  if  we  assume  symmetric  pulse  shapes,  then  it 
can  be  shown  that  [9] 


I 


+1 


(3.31) 


and 


P(e)  -  }!PAI.H0-O,d.1.l,dtl.ll 

+  T6(i>c<e|V0'd-i=0'd+i=1) 

+  5ipGlelV°'1-r°'d*f01 


+  PB(aiV1'd-1=l'd-H*l) 

+  PD 1 d0*d ' d-l*P' d+l=1^  ^ 

+  PH(e|dQ=l,d_1=0,a+1=0) ]  . 


(3.32) 


3 . 5  Gaussian  Approximation  for  Shot  Noise  Distribution 

If  we  approximate  the  shot  noise  distribution  as  Gaus¬ 
sian,  the  mean  and  variance  of  the  process  completely  des¬ 
cribe  it.  Further,  the  difficulty  in  evaluating  an 
inhomogeneous  nonstationary  shot  noise  process  is  avoided. 

If  we  assume  the  avalanche  gain  of  the  APD  is  large,  the 
Gaussian  distribution  is  a  good  approximation  to  the  exact 
distribution  [10] .  It  can  be  shown  [9]  that  the  expression 
for  the  average  probability  of  error,  P(e),  in  (3.32)  can  be 


expressed  as 


p<e»  -  5  +  16 


erf 


/  a  -  ( 1+2AI ) 

- t -  SNR  -  erf 


aT  -  2AI 


\  /2<1  + 


2AI , 
N+l 


'2(1  + 


2AI-1' 

N+l 


♦* 


erf 


am  -  ( 1+AI ) 


AT 

/2(1  +  STI 


SNR  -  erf 


aT  -  AI 

N+l  ; 


SNR 


16 


erf 


/  a„  -  1 

1 

T* 

SNR 

-  erf 

T 

SNR 

,<2N_ 

1 

/ 

/ 

V  /  N+l 

where  dimensionless  quantities  have  been  introduced.  A 
measure  of  the  ISI  is  given  by  the  spreading  ratio 


SNR 


(3.33) 


AI  = 


(3.34) 


We  define  the  signal  energy  by 


s  =  rr  h  I 
hf  o 


(3.35) 


The  decision  threshold,  normalized  by  the  signal  energy,  is 
given  by 

rT 

<*T  =  i1  •  (3.36) 

The  ratio  of  thermal  noise  to  shot  noise  is  given  by 

thermal  noise  variance 


N  = 


ath2 


hs 


shot  noise  variance 


(3.37) 
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We  define  the  signal- to-noise  ratio  (SNR)  as 
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CHAPTER  4 

RECEIVER  OPTIMIZATION 


4 . 1  Evaluation  of  the  Optimum  Decision  Threshold 

If  the  receiver  is  operating  under  the  ideal  conditions 
of  zero  thermal  noise,  zero  ISI,  and  zero  dark  current,  the 
resulting  decision  threshold  would  be  equal  to  zero.  A 
given  time  slot  will  contain  energy  only  if  a  one  is  trans¬ 
mitted  in  that  period.  There  will  be  no  energy  present  in 
the  time  slot  if  a  0  is  transmitted.  Unfortunately,  both 
thermal  noise  and  ISI  are  present,  forcing  the  decision 
threshold  to  be  greater  than  zero.  By  optimizing  this  deci¬ 
sion  threshold,  we  can  improve  on  error  performance  of  the 
receiver . 

A  measure  of  the  ISI  is  the  spreading  ratio,  II,  dis¬ 
cussed  in  the  previous  chapter.  If  we  assume  Gaussian  pulse 
shapes,  expressions  can  be  easily  derived  for  the  pulse  ener¬ 
gies  IQ  and  1^  in  terms  of  the  error  function  and  complemen¬ 
tary  error  function  [9] .  The  analytical  expressions  for  the 
pulse  energies  are 


I  =  erf 
o 


2/2  T 


out 


(4.1) 


1^  =  j  erfc 


2/2  T 


(4.2) 


out  / 


1 


I 

) 


We  can  express  the  spreading  ratio  in  terms  of  Ts 

by 


AI 


~  erfc 


2  /2  T 


out 


erf 

1  Ts  \ 

2 ^  T  „ 

V  °ut  j 

and  T 


out 


(4.3) 


In  Figure  8,  the  spreading  ratio  is  plotted  as  a  function  of 
the  ratio  of  signaling  interval,  Ts,  to  the  rms  output  pulse 
width,  T  ,  for  a  Gaussian  pulse  shape. 

To  determine  the  optimum  threshold,  we  will  examine  the 
two  conditions  under  which  a  detection  error  is  most  likely 
to  occur.  The  receiver  is  most  likely  to  erroneously  detect 
a  1  in  Tsq  when  both  adjacent  time  slots  contain  l's.  Energy 
will  spread  into  Tsq  from  adjacent  time  slots  due  to  ISI. 
Thermal  noise  will  also  increase  the  energy  in  Tsq.  If  the 
ISI  and  thermal  noise  contributions  are  great  enough,  the 
decision  threshold  will  be  exceeded  and  the  receiver  will 
erroneously  detect  the  presence  of  a  1  in  Tsq. 

The  receiver  is  most  likely  to  decide  a  0  was  trans¬ 
mitted  in  Tsq  when  a  1  was  actually  transmitted  when  the 
adjacent  time  slots  both  contain  0's.  If  the  ISI  is  severe 
enough,  excessive  energy  will  spread  into  Ts_^  and  Ts+^, 
leaving  an  insufficient  amount  of  energy  in  TsQ  to  exceed 


the  threshold. 


Figure  8.  ISI  spreading  ratio  as  a  function  of  the  ratio  of  the 
signaling  interval  T  to  the  rms  output  pulse  width  T 
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An  approximation  for  the  probability  of  error  can  be 
found  if  these  two  worst  case  conditions  are  summed  together 
in  the  expression 


P(e)  *  PA(e|do  =  0,  d_x  =  1,  d+1  =  1) 

+  PH(e|do  =  d-l  =  °'  d+l  =  0)- 


(4.4) 


To  determine  the  optimum  threshold  which  minimizes  P(e),  we 
will  differentiate  (4.4)  with  respect  to  the  normalized 
threshold,  a^,  set  the  result  equal  to  zero,  and  solve  for 
the  optimum  threshold,  aQpT*  The  expression  for  the  optimum 
threshold  is 

_  N ( 2  -  4AI)  -  (ab  +  c  +  d  +  e) 1/2  ,  ,  c, 

^OPT  ~  4AI  -  2  1 

where 

a  =  AI(16AI  +  8N  -  8)  -  4N 

2AI  -  1  2 

b  =  ln(l  +  )  /  SNI< 

c  =  N2(16AI2  -  16AI  +  4) 
d  =  N ( 3 2  AI3  -  16AI2  -  8AI  +  4) 
e  =  32AI 3  -  32AI2  +  8AI  . 

We  see  that  in  the  absence  of  ISI  and  thermal  noise,  the 
optimum  threshold  reduces  to  zero  as  expected. 
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4 . 2  Numerical  Results 

In  Figures  9-12,  the  optimum  threshold,  aopT,  is  plotted 
as  a  function  of  the  spreading  ratio,  LI,  for  several  values 
of  the  thermal  noise  to  shot  noise  ratio,  N.  Values  of  N 
used  are  0.0001,  1.0,  2.0,  and  3.0  in  each  figure.  The  SNR 
is  held  constant  for  each  set  of  curves.  The  advantage  of 
increasing  the  signal  strength  is  evident  if  we  compare 
Figure  12  with  Figure  9  and  note  the  significantly  lower 
values  of  optimum  threshold  when  a  strong  signal  is  present. 

In  Figures  13-16,  the  optimum  threshold  is  plotted  as 
a  function  of  the  spreading  ratio  for  several  values  of  SNR. 
The  value  of  the  thermal  noise  to  shot  noise  ratio,  N,  is 
held  constant  for  each  set  of  curves.  For  values  of  N 
greater  than  0.5,  the  value  of  the  optimum  threshold  in¬ 
creases  linearly  with  increasing  values  of  LI.  For  small 
values  of  N  the  relationship  between  '*QpT  and  Al  is  very 
nonlinear.  Under  these  conditions,  very  small  optimum 
threshold  values  can  be  realized  for  low  ISI  levels. 


Figure  10.  Optimum  threshold  as  a  function  of  the  ISI  spreading  ratio 
for  a  Gaussian  shot  noise  distribution.  SNR  is  held 


Figure  11.  Optimum  threshold  as  a  function  of  the  1SI  spreading  ratio 
for  a  Gaussian  shot  noise  distribution.  SNR  is  held 


Figure  12.  Optimum  threshold  as  a  function  of  the  ISI  spreading  ratio 
for  a  Gaussian  shot  noise  distribution.  SNR  is  held 


Figure  13.  Optimum  threshold  as  a  function  of  the  IS1  spreading  ratio 
for  a  Gaussian  shot  noise  distribution.  N  is  held  constant 


Figure  14.  Optimum  threshold  as  a  function  of  the  ISI  spreading  ratio 
for  a  Gaussian  shot  noise  distribution.  N  is  held  constant 
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Figure  15.  Optimum  threshold  as  a  function  of  the  ISI  spreading  ratio 
for  a  Gaussian  shot  noise  distribution.  N  is  held  constant 


Figure  16.  Optimum  threshold  as  a  function  of  the  1SI  spreading  ratio 
for  a  Gaussian  shot  noise  distribution.  N  is  held  constant 


CHAPTER  5 


RECEIVER  PERFORMANCE  AND  SYSTEM  DESIGN  CONSIDERATIONS 

5 . 1  Receiver  Error  Rate  Performance 

Since  we  have  expressions  for  the  probability  error 
(3.33),  and  the  optimum  decision  threshold,  (4.5),  we  an 
replace  aT  in  (3.33)  by  ^0PT  to  obtain  an  expression  for 
P(e)  which  is  minimized  for  the  variables,  AI,  SNR,  and  N. 

If  the  communications  system  is  designed  such  that  the  spread¬ 
ing  ratio,  AI,  is  zero,  P(e)  becomes  a  function  of  N  and  SNR 
only  and  is  shown  in  Figure  17.  The  error  probability  is  a 
monotonically  decreasing  function  of  the  signal- to-noise 
ratio.  For  a  constant  signal- to-noise  ratio,  the  error 
probability  increases  with  increasing  values  of  N  which  cor¬ 
responds  to  an  increase  in  the  thermal  noise.  Also,  to  main- 

_  g 

tain  a  specified  error  rate  of  10  ,  the  SNR  must  increase 

by  approximately  2.5  aB  if  the  noise  ratio  increases  from 
lO-4  to  10.0. 

In  Figures  18-21,  P(e)  is  plotted  versus  AI,  SNR,  and 

N  when  the  intersymbol  interference  is  not  zero.  Again,  the 

error  probability  decreases  monotonically  as  a  function  of 

the  signal- to-noise  ratio.  N  is  held  constant  for  each  plot 

-4 

and  takes  on  values  of  10  ,  0.1,  0.5,  and  10.0  for  succes¬ 

sive  figures.  In  each  plot,  P(e)  is  shown  for  I  =  0.01, 

0.1,  0.2,  and  0.3.  A  review  of  Figures  18-21  indicates  that 


± 


Average  probability  of  error  as  a  function  of  the  signal 
to-noise  ratio  for  various  values  of  the  spreading  ratio 
with  N  held  constant. 
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the  spreading  ratio  has  a  much  greater  effect  upon  the  error 
probability  at  low  noise  levels.  Thus,  when  designing  an 
optical  fiber  communications  system  with  a  specified  error 
rate,  controlling  the  ISI  becomes  more  important  at  low 
noise  levels. 

5 . 2  System  Design  Considerations 

When  designing  an  optical  fiber  communications  system, 
the  error  probability  is  of  primary  importance  and  is  normally 
one  of  the  initial  specifications.  A  typical  value  of  N  can 
be  obtained  for  available  receivers.  Knowing  the  error 
probability  and  N,  and  by  referring  to  Figures  18-21,  a  maxi¬ 
mum  value  for  AI  can  be  obtained  once  we  specify  an  appro¬ 
priate  SNR  based  upon  the  system  parameters.  In  Figure  8, 

we  can  obtain  a  value  for  T  /T  .  for  this  value  of  AI. 

s  out 

Specifying  the  bit  rate  now  allows  us  to  calculate  the  value 

for  T  . .  Knowing  the  maximum  allowable  value  of  T  and 

out  out 

the  desired  optical  fiber  cable  length  further  specifies 
the  value  for  the  rms  input  pulse  width.  We  can  now  select 
the  proper  fiber  and  source  in  order  to  meet  the  requirements 
for  the  optical  fiber  system. 
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CHAPTER  6 
CONCLUSION 


In  this  report  we  have  analyzed  the  effects  of  pulse 
dispersion  and  intersymhol  interference  on  the  performance 
of  a  graded-index  fiber  optic  communication  system.  By 
using  the  WKB  approximation,  it  was  found  that  with  the  proper 
choice  of  the  index  gradient,  a,  modal  dispersion  can  be  re¬ 
duced  such  that  material  dispersion  is  the  dominant  factor 
in  calculating  pulse  dispersion  in  a  graded-index  fiber. 

An  expression  for  the  rms  output  pulse  width,  TQut, 
was  developed  in  terms  of  the  dispersion  parameters,  the 
length  of  the  fiber,  and  the  rms  input  pulse  width,  T^.  It 
was  shown  that  a  narrow  linewidth  source,  such  as  a  laser, 
is  capable  of  achieving  a  smaller  TQut  than  an  LED  with  an 
equivalent  value  for  .  We  also  find  that  for  small  values 
of  T^,  the  pulse  underwent  severe  temporal  spreading,  which 
resulted  in  TQut  being  many  times  greater  than  . 

A  mathematical  model  of  an  optical  digital  receiver 
was  developed  using  an  APD  as  the  optical  detector  and  an 
integrate-and-dump  pulse  detection  scheme.  Dark  current  was 
considered  negligible  and  a  Gaussian  approximation  was  made 
for  the  shot  noise.  An  expression  for  the  optimum  decision 
threshold  of  the  receiver  was  developed  in  terms  of  the  SNR, 
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the  thermal  noise  to  shot  noise  ratio,  and  the  spreading 
ratio,  AI.  We  have  shown  that  a  linear  relationship  exists 
between  aQp T  and  AI  when  N  is  greater  than  1.0.  In  general 
we  have  shown  that  c*0pT  increases  with  increasing  values  of 
N  and  AI.  Also,  ot^p^  decreases  with  increasing  values  of 
SNR. 

An  expression  for  the  average  probability  of  error 
was  developed  for  the  system  as  a  function  of  the  optimum 
decision  threshold,  the  spreading  ratio,  the  thermal  noise 
to  shot  noise  ratio,  and  the  SNR.  In  general  the  error  rate 
increases  for  increasing  values  of  N  and  ISI,  while  decreas¬ 
ing  for  increasing  values  of  the  SNR.  Also,  ISI  has  a  more 
pronounced  effect  on  the  error  rate  at  low  noise  levels 
when  compared  to  the  situation  at  high  noise  levels. 

In  the  final  section,  we  developed  a  method  of  design¬ 
ing  a  fiber  optic  communication  system  based  upon  the  pre¬ 
ceding  analysis.  Using  this  method,  an  appropriate  selec¬ 
tion  of  the  optical  source  and  fiber  optic  cable  can  be  made 
based  upon  the  system  specifications. 
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